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Abstract Wc present a new approach to combine remote observations and in- 
situ data by STEREO/HI and Wind, respectively, to derive the kinematics 
and propagation directions of interplanetary coronal mass ejections (ICMEs). 
We use two methods, Fixed-</> (F</>) and Harmonic Mean (HM), to convert 
ICME elongations into distance, and constrain the ICME direction such that 
the ICME distance-time and velocity-time profiles are most consistent with in- 
situ measurements of the arrival time and velocity. The derived velocity-time 
c"| functions from the Sun to 1 AU for the three events under study (1-6 June 

2008, 13-18 February 2009, 3-5 April 2010) do not show strong differences 
for the two extreme geometrical assumptions of a wide ICME with a circular 
front (HM) or an ICME of small spatial extent in the ecliptic (Fcf>). Due to the 
geometrical assumptions, HM delivers the propagation direction further away 
from the observing spacecraft with a mean difference of ~ 25°. 



1. Introduction 

Coronal mass ejections (CMEs) are manifestations of the most powerful erup- 
tions on the Sun and are expulsions of a huge amount of plasma and the 
embedded magnetic field. Their velocities range between a few 100 km s _1 up to 
more than 3000 km s _1 . The frequency of occurrence correlates with the solar 
cycle, and faster and more powerful events are more common during the solar 
maximum phase. They can have wide latitudinal and long itudinal exten t s and 
their typical masses are of the order of w 10 12 - 10 13 kg (c/. iGosling et ail . \l974 
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Hundhausenl . 1 1 9971 ) . CMEs propagating from the Sun through the interplanetary 



space are called interplanetary CMEs (ICMEs). 

Before 2006 it was not possible to directly link CMEs and their properties as 
measured in-situ. A milestone in the investiga tion of CMEs is the NASA Solar 
TErrestrial RElations Observatory (STEREO: iKaiser et al. . 20081 ) mission with 
its twin satellites STEREO- AHEAD (A) and STEREO-BEHIND (B). STEREO- 
A is leading the Earth in its orbit around t he Sun and STERE O-B is following. 
The Heliospheric Imagers (HI1 and HI2: Evles et al. . 2009J) arc part of the 
Sun-Earth Connection Corona l and Heliospheric Inv estigation instrument suite 
fSECCHI: lHoward et all 120081: lHarrison et aZ.l . l2009l ) onboard the two STEREO 
spacecraft and they enable us for the first time to observe solar transient events 
from two vantage points from outside the Sun-Earth line. Such observations 
constitute a unique way of investigating the behavior of ICMEs all the way from 
the Sun to Earth. In addition, the possibility to study solar minimum events is a 
big advantage because it is necessary to do case studies of CMEs showing simple 
and well-defined remote as well as in-situ signatures. 

There arc different methods that can be used to infer the direction of prop- 
agation of an ICME. So me methods use single spacecraft ob s ervat ions, e.g. the 
Fixed-0 fitting method dSheelev efall Il999l: iRouillard et~aZI. l2008h or the Har- 
monic Mean fitting procedure ( Lugad . 2010l : Mostl et all . 201 lfl while o t hers us e 
data from both spacecraft, e. g. the triangulation method bv lLiu et al. ( 2010a ). 
extended to circular fronts bv lLugaz et all (|2010l ). 

The aim of this work is to analyze the kinematics and propagation directions of 
a set of CMEs up to a distance of 1 AU. For this we developed a method based on 
a combination of remote and in-situ measurements. We constrain measurements 
from time-elongation plots (Jmaps), produced out of the white-light HI1/HI2 
images along the ecliptic plane, with the in-situ measured arrival time and arrival 
velocity The best matches deliver propagation directions and velocity profiles 
based on the geometric assumptions we use, for either very wide or narrow ICME 
fronts. This should serve as a basis to develop adequate methods to predict the 
arrival times of Earth-directed CMEs. 



2. Data 



The Heliospheric Imagers (HI: lEvles et aZ.l . l2009l ) onboard STEREO for the first 
time give us the possibility to perform remote sensing in white-light between 
the Sun and the Earth from outside the Sun-Earth line. Because of their wide 
observation angles (HI1: 4-24°, HI2: 18.7 - 88. 7°) and scattering effec t s, the 



interpretation of these images is challenging (e.g. Vourlidas and Howard . 20061 : 



Kahler and Webbl . l2007t iHoward and Tappinl . l2009h . To derive the kinematics of 
ICME fronts on their way through the inner heliosphere we used remote observa- 
tions of both HI instruments on STEREO-A as well as the in-situ p roton density 



19951 ) and STEREO-B (PLASTIC: [Galvin et aZ.l . l2008l ) 



pr 

and velocity delivered by the Wind spacecraft ne ar Earth (SWE: iQgilvie et al 



In this study, we discuss three CME-ICME events observed end-to-end from 
the Sun to 1 AU, covering a wide range of CME initial conditions and param- 
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eters. The first event of 2-6 June 2008 is an example of a v ery slow streamer- 
blowout type CME with a small longitudinal angular width dMostl et al. L 2009; 



Robbrecht. Patsourakos. and VourlidasLl200"i Lynch et aZ.l . l2010HWood. Howard, and Sock"erl | 
2010h . The CME of 13 - 18 February 200 9 originated from a bipolar active region 
and was associated with an EUV wave ( Kienreich. Temmer. and Veronid . 20091 : 
Patsourakos and Vourlidasl 120091 ). It was slow as well b ut reached its propaga - 



tion velocity of about 350 km s very close to the Sun (IMiklenic et aZ.1 . 1201 lh . 



The event is ex pected to have a w ider extent in the ecliptic because of a lower 
axis inclination (jMostl et all 120111 ). The third eve nt of 3-5 April 2010 was the 



first f ast and geoeffective ICME of Solar Cycle 24 (jMostl et aill2010l : iLiu et al 



20111 : IWood et all . 1201 lh and led to a damage of the Galaxy 15 satellite in 



geosynchronous orbit and to what has been called a perfect substorm. 



3. Methods 

For the interpretation of HI images it is necessary to consider scattering effects. 
ICMEs are detected as the photospheric light s cattered off free electrons in 
the C ME body, called Thomson scattering (e.g. Billings . 19661 : iHundhausenl . 



19931 ). The scattered light has its maximum when the line of sight (LoS) is 
perpendicular to the line between the Sun and the scattering particle. That is 
the reason light that we see in the white-light images originates on a sphere 
with the Sun-observer l ine as its diameter, also called the Thomson-surface (TS: 
Vourlidas and HowardLl2006h. W hen le aving the TS the intensity of the scattered 



light decreases (e.g. iMorrill et aLl . l2009h . Brightness dis tribution and morphology 



arc different depending on the ICME axial orientation (jCremades and Bothmer 



20041 ) . Thomson scattering influences the geometrical derivations of the following 



techniques. 



Sig natures of ICMEs ar e measured within time-elongation plots fJmaps: ISheelev et oi.l .| 



19991: iDavies et all l2009h which are produced out of stripes along the ecliptic 



plane extracted from heliospheric images. These stripes are rotated and aligned 
next to each other with the time increasing in the x- and the elongation in the 
y-dircction. When an ICME feature is measured within a Jmap the distance 
from Sun-centre is denoted in elongation, i.e. in angular degrees. To use the 
measurements for further analysis we arc interested in the ICME kinematics 
expressed in radial distance, and thus the measured elongation has to be con- 
verted. These calculations are limited because of the supposed shape of the 
ICME front, which cannot be determined by using remote sensing from only 
on e vantage point. These methods and techniques are reviewed and discussed 
in Liu et al. ( 2010bh . In the following calculations we converted the elongation 



angle into radial distance by using two methods that make different geometrical 
assumptions for the shape of the ICME front: Fixed-</> and Harmonic Mean. 

3.1. Fixed-0 Method 

The simplest way to c onvert elongation into distance is the Point P (PP) method 
( Howard et al. . 20061) . It assumes a CME as a circle all the way around the Sun. 
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Figure 1. Comparison between Point P, Fixed-</>, and Harmonic Mean methods. The gray 
line shows the line of sight, angle e is the measured elongation, d is the Sun-observer distance, 
(j> is the propagation direction relative to d . For the same value of <j>, Fixed-0 gives a larger 
distance than HM. 



In co ntrast to PP the Fixcd-0 method (F6: ISheelev et all 1 1 9991 : iKahler and Webbl .| 



20071 ) assumes a radial propagation of a single plasma element along a straight 
line (see Figure [TJ, i.e. a constant propagation direction. Thi s is an important 



difference to other methods (e.g. the triangulation technique bv lLiu et qLl . l2010al) 



that make geometrical assumptions of a point wise extent of the CME as well but 
use two different vantage points (STEREO-A and STEREO-B) and are therefore 
able to determine the direction for each point along the track. This approach 
converts the measured elongation angle into radial distance from the Sun for a 
given propagation direction of the ICME: 



RF(j>(t) 



sine(i) 
°sin(e(i) + </>) : 



(1) 



where Rp^it) is the calculated distance, d is the Sun-observer distance, e(t) the 
measured elongation angle and (j) the derived propagation direction, measured 
away from the observer, with positive values corresponding to solar West. 

3.2. Harmonic Mean Method 



For wide CMEs, the Harmonic Mean method (HM: iHoward and Tappinl . 120091 ; 
Lugaz. Vourlidas. and Roussev , 2009) may be more appropriate than F</>. It as- 
sumes that the measured part is not a single particle but a sphere (or a circle in 
the ecliptic plane) connected to the Sun at all times. As shown in Figure \T\ it 
further assumes that the observer always looks along the tangent to this circle. 
The resulting equation can be understood as the harmonic mean of the PP and 
F(f> methods: 
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Figure 2. Resulting radial distance as function of time for five different propagation directions 
for the F<f> (upper panel) and the HM (lower panel) methods. The angle is the propagation 
direction relative to the observer — negative means eastward. The solid horizontal lino marks 
the location of the spacecraft providing the in-situ measurements and the vertical solid line 
indicates the arrival time of the ICME at the in-situ spacecraft. 



^ = T^7& " (2) 

1 + sin(e(t) + </>) 

where RiiM{t) is the calculated distance, d a is the Sun-observer distance, e(i) 
the measured elongation angle, and <j> the derived propagation direction. This 
was introduced bv lLueaz. Vourlidas. and Roussev ( 2009t) who found that CME 



velocities at large elongation angles in a simulation were reproduced better with 
the HM method, and were over- and underestimated by the Fixed-$ and Point P 
methods, respectively. 
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3.3. Combining Remote and In-Situ Signals 

We have developed a new technique, based on the two conversion methods, by 
adding two boundary conditions: first, the arrival time [t a ] of the ICME front 
at the in-situ spacecraft and second, the in-situ measured velocity [V^]. As an 
illustration, Figure [5] shows the resulting distance-time profiles of the 1-6 June 
2008 event calculated with different propagation angles (Top: F</>, Bottom: HM). 
The cross is the measured arrival time of the ICME front at the location of the 
in-situ spacecraft. 

3.3.1. Boundary Condition: In-Situ Arrival Time 

For estimating the most reliable propagation angle [</>] the measured elongation 
value at the in-situ arrival time of the ICME [e t J is converted with F</> into 
a distance, RF<p,t a {4>), by using different angles (cp G [0°,180°]). i?F0,t a (0) is 
subtracted from the distance between the Sun and the in-situ spacecraft [di]: 

Ad F(t ,(^) = Rp^tM) ~ di- (3) 

The angle which minimizes A.dF<f> (</>), which we call <fidF(j>, is the resulting direc- 
tion from the constraint with the arrival time. 

To do the same calculation for HM we have to consider the particular, circular 
geometry of the front. To this end one cannot use the whole diameter of the circle 
[Rhm], but rather the distance between the Sun and the point of intersection 
of the HM circle with the line connecting the Sun and the in-situ spacecraft, as 
illustrated in Figure [3] We call this distance RiiMi,t a (4>), and it is calculated by 

RHMLt a O) = RHM,t a (0) cos S, (4) 
where <5 is the angle between di and Rhm- Similar to above, we calculate: 

Ad H Mi(4>) = RnMi,t a {4>) ~ di, (5) 

where AdHMi(4>) 1S the resulting difference and RHMi.t a {4>) the result of HM 
along di. The propagation direction [4>dHAi] corresponding to the minimum of 
AdHMi(4>), is again the outcome of the technique. 

3.3.2. Boundary Condition: Velocity at In-Situ Arrival Time 

We applied the same minimization method as described before by subtracting 
the in-situ measured velocity [Vi] from the calculated velocity converted from 
the measured elongation value at t a , VF0,t a (0 £ [0°, 180°]): 

&V F<t> {$) = V F<t> ,tM-Vh (6) 

where AVf4,(4>) is the calculated difference, VF4,.t a {4>) the derived velocity at 
arrival time, and Vi is the in-situ measured velocity of the ICME front. 
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Figure 3. Illustration showing the calculation of the distance [Ad] from the Sun to the point 
of the circular CME front in the direction of the in-situ observing spacecraft. The gray line 
shows the line of sight, angle e is the measured elongation, d is the Sun-observer distance, <fi 
is the propagation direction relative to d a , Rhm is the distance of the CME apex from the 
Sun, RjjMi is the distance of the CME front in direction of the in-situ spacecraft from the 
Sun, and 8 is the angle between Rhm and Rhah- 



The same method can be applied by using HM. The velocity Vhm% was 
calculated on the basis of Rhah, thus in the direction of the in-situ spacecraft, 
which for HM does not necessarily have to be the apex of the ICME front, but 
can, in principle, be also any point along the circle. We define 



AVhah (0) = V H Mi,u V u (7) 

where AVHMi(<f>) is the calculated difference and VjiMi,t a (</>) the derived velocity 
at arrival time. The angle belonging to the minimum difference [4>vf<s> or 4>vhm] 
is the dir e ction r esulting from the constraint with the velocity at arrival time. 
Liu et all (|2010bl ) pointed out that there is an ambiguity for possible propagation 



directions for HM. Because of mathematical reasons there are two minima within 
the range of £ [0°, 180°]. Since we get two minima for both constraints (Ad, 
AV) we choose the value where both minima show the best agreement with 
each other; in our cases, it seems the smaller value. The conversion using the 
next minimum would yield a large distance of the CME apex, which may not be 
possible given the observed CME size from two vantage points (Figure 0]) . 

The directions from the arrival time and velocity constraints will not be 
identical. To combine the resulting propagation directions for each method, we 
define 
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Figure 4. Two possible solutions for HM. The green arrow shows the direction of the ICME 
apex resulting from the first minimum and the green circle shows the ICME front (CME 1) for 
this direction at arrival time. The red arrow indicates the direction derived from the second 
minimum of AV for the ICME event of June 2008. The apex of this ICME (CME 2) would be 
2.5 AU away from the Sun at arrival time at the in-situ spacecraft. 



4>F<f> = {(pdFij, + <t>VF<f>)/2 (8) 
<pHM = (fidHM + <?VffM)/2 (9) 

so we take the average values for the direction for each method and use these 
for further analysis. Here, 4>hm is the propagation direction of the ICME apex 
for the HM method. However, all kinematics are calculated for the part of the 
circle which hit the in-situ spacecraft and not for the apex. 

The distance-time profile for the defined mean propagation angle was fitted 
by using a cubic spline, from which we derived the velocity profile via numerical 
differentiation. All errors result from the manual tracking of the ICME front only. 
To be able to do an error estimation, every feature was measured five times. The 
resulting standard deviation also yields a reliable error for the numerical velocity 
derivation using three-point, Lagrangian interpolation. We also indicate an error 
for the propagation direction for both methods. 
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4. Results 



Three different ICMEs were investigated to cover a range from slow to fast CME 
initial speeds. One event was slow and narrow (June 2008), the second was also 
slow but had a wide longitudinal extent (February 2009), and the third ICME 
(April 2010) was a fast and gcocffcctivc event. 

4.1. 02-06 June 2008 ICME 

The leading edge of this ICME was seen on 2 June 2008 09:00 UT in the HI1 
field of view (FoV) and on 3 June 2008 12:00 UT in the HI2 FoV of STEREO- 
A. This CME had no obvious signatures of m agnetic reconnection on the Sun 



( Rob brecht. Patsourakos. and VourlidasL 120091) . A magnetic cloud was detected 



between 6 June 2008 sa 22 UT and 7 June 2008 » 12 UT at STER EO-B 
(jMostl et all 120091 iLvnch et all |201CA IWood. Howard, and Sockerl . |201dh . The 



separation between both STEREO spacecraft was 53.9° at that time. Figure [5] 
shows the proton density (STEREO-B, PLASTIC) and the Jmap (STEREO-A, 
HI1/2) with two clear traces of the event. The first one is the leading edge and 
the second one is the core of the ICME. Both Jmap-tracks match with the in-situ 
arrival times (sudden increase of the proton density at STEREO-B) as indicated 
with the vertical red lines. Its in-situ signature as well as the heliospheric images 
both reveal a pronounced three-pa rt structure, with the density enhanc ements 
bracketing the magnetic flux rope (jMostl et all 120091: ILvnch et all |2010[ ). 



To obtain the propagation direction and the velocity only the first track, i.e. 
the ICME leading edge, was investigated. Figure [5] shows the calculation of the 
minima of the differences between the converted distance at in-situ arrival time 
and the distance between Sun and in-situ spacecraft (4>dF<j> and (pdHu) as well 
as for the differences between the derived velocity at arrival time and the in-situ 
measured arrival velocity [4>vf<i> and 4>vhm}- ^4> (top panel) yields as a mean 
direction E24±3° (relative to the Sun-Earth line) and HM (bottom panel) results 
in E51 ± 6°. 

Figure [7] shows the kinematics of this event derived from F(f> using the ob- 
tained propagation angle. The upper panel displays the result of the conversion 
from elongation into solar radii and the link to the arrival time at the in-situ 
spacecraft. The standard deviation lies between 0.5 R© (HI1) and 1.6 R© (HI2). 
The middle panel shows the direct derivation of the measurements (crosses) and 
the derivative of the spline fit. The kinematics using F(j> (4>F4> = E24 ± 3°) yield 
a continuous acceleration from w 330 to m 440 km s _1 and a mean velocity 
of 389 ± 48 km s" 1 . The velocity at arrival at STEREO-B derived from Fcj> 
is 432 ± 61 km s _1 . F<p results in a typical late acceleration at w 1 AU. This 
can not be interpreted as a tru e behavior of the CME but is an arte fact of the 
method (cf. the d iscussions in Lueaz. Vourlidas. and Roussev . 20091 . and) and 



Wood et al\ |2009) . The new method using in-situ data as constraints diminishes 
this effect. The same approach but for HM is shown in Figure |U The errors of 
the measurements are between 0.4 R Q (HI1) and 1.2 R Q (HI2). Kinematics 
using HM (4>hm = E51 ± 6°) deliver a mean velocity of 395 ± 35 km s _1 . The 
CME shows a strong acceleration from 280 to « 440 km s _1 up to a distance 
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Figure 5. In-situ measured proton density from STEREO-B (top). The first red line from 
the left marks the time of the ICME arriving at STEREO-B and others delimit two strong 
peaks in the proton density. The lower panel shows the Jmap produced from remote sensing 
data of hcliosphcric images of STEREO-A with overplottcd measurement points (red crosses) . 
The white horizontal line indicates the position of STEREO-B. 



of about 100 solar radii followed by a slight deceleration to a final velocity of 
422 ± 44 km s _1 , consistent with the in-situ measured velocity at STEREO-B 
of 403 km s" 1 . 



4.2. 13-18 February 2009 ICME 



The ICME arrived on 13 February 2009 11:00 UT in the HI1 FoV of STEREO- 
A and became visible in HI2 on 14 Februar y 2009 04:00 UT. This CME was 
assoc iated with a flare and an EUV wave ( Kienreich. Temmer. and VeronigL 
2009) occurring at the limb as seen by STEREO-A, since the separation of the 
STEREO satellites was 91° at that time. Compared to the June 2008 event there 
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Figure 6. Propagation directions for Fixed-0 (top panel) and Harmonic Mean (bottom panel). 
The black crosses show the differences between the calculated distance at arrival time and the 
distance of the in- situ spacecraft from Sun-centre for different propagation angles. The red 
asterisks indicate the same approach but for the difference of the calculated velocity and the 
in-situ measured velocity at arrival time. 



is no characteristic three-part structure — neither in the heliospheric images nor 
in the in-situ data. STEREO-B measured an incr ease of the p r oton density that 
was followed by a large-scale magnetic flux rope ( Mostl et al . 201 lh . 



Figure [5] shows the proton density measured by STEREO-B and the Jmap 
from STERE O-A. The ICME trac k became faint rather quickly, as expected for 
a limb CME ( Morrill et all . l2009h . Comparing the track in the Jmap and the 
proton density in Figure shows that the ICME seems to arrive earlier at the 
elongation of STEREO-B than it was measured in-situ. 

Figure [TTJ] illustrates the calculation of the propagation directions = 
E36 ± 4°. 4>hm = E61 ± 3°). Again. HM gives a direction further away from the 
observer. The upper panel in Figure ITTI shows the result of the conversion from 
elongation into radial distance [F<^>] and the linkage with the in-situ arrival time. 
The standard deviation of both conversion methods is between 0.5 R Q (HI1) 
and 1 Rq (HI2). Kinematics using F<f> yield a mean velocity of 310 ± 29 km s _1 . 
The velocity profile shows a nearly constant speed up to a distance of » 100 
Rq followed by an acceleration from w 280 km s _1 to a velocity at arrival time 
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Figure 7. Top: Resulting distance-time profile for F<fi. The crosses are the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal 
line shows the arrival time of the ICME at STEREO-B. Middle: The solid curve shows the 
derivation of the fit, which is done to determine the velocities of this event. The gray area 
indicates the standard deviation of the measurements. The horizontal line indicates the in-situ 
measured velocity at STEREO-B. Bottom: Residuals of the fit and the direct measurements. 



of 416 ± 35 km s -1 . Again F</> shows an apparent acceleration close to 1 AU. 
The mean velocity derived of HM is 325 ± 28 km s _1 and the impact velocity 
is 398 ± 34 km s _1 , whereas the in-situ measured impact velocity of the ICME 
at STEREO-B is ~ 362 km s -1 (sec Figure IT2"j) . There seems to be an upward 
kink in the track in the Jmap at about 35° elongation that is also visible in the 
derived kinematics as an acceleration. 

4.3. 03-05 April 2010 ICME 

Compared to the previous mentioned events this CME was relatively fast for 
solar minimum. The ICME arrives in the HI1 FoV of STEREO-A on 3 April 2010 
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Figure 8. Top: Resulting distance-time profile for HM. The crosses are the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal 
line shows the arrival time of the ICME at STEREO-B. Middle: The solid curve shows the 
derivation of the fit, which is done to determine the velocities of this event. The gray area 
indicates the standard deviation of the measurements. The horizontal line indicates the in-situ 
measured velocity at STEREO-B. Bottom: Residuals of the fit and the direct measurements. 



10:00 UT and a shock was detected on 5 April 2010 07:5 8 UT at Wind, situated at 



2010h . The separation between 



the LI point sunward of the Earth (IMostl et al. 
STEREO-A and Wind was 67.4° at that time. It was a geoeffective ICME and 
caused a moderate geomagnetic storm with a minimum of the Dst ind ex about 
-72 nT and maximum Kp of 8 (jMostl et all l201(A IWood et aZ.l . l201lh . 

Figure [T3] shows the proton density at Wind and the Jmap of HI1A and 
HI2A of the event. The track of the ICME fits well with the sharp increase in 
proton density. As illustrated in Figure [JJ] the two conversion methods yield a 
propagation direction of = W3 ± 4° and 4>hm = E25 ± 10°, respectively. 
The measurement errors for F(j> are 1.3 R© (HI1) and 1.5 R© (HI2). for HM 
1 R© (HI1) and 1.3 R© (HI2). Kinematics from F(f> result in a mean velocity 
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Figure 9. In-situ measured proton density from STEREO-B (top). The first red line from 
the left marks the time of the ICME arriving at STEREO-B and others delimit two strong 
peaks in the proton density. The lower panel shows the Jmap produced from remote sensing 
data of hcliospheric images of STEREO-A with overplotted measurement points (red crosses) . 
The white horizontal line indicates the position of STEREO-B. 



of 829 ± 122 km s" 1 and an impact velocity of 825 ± 129 km s" 1 (Figure ITS) . 
kinematics from HM yield a mean velocity of 854 ± 100 km s _1 and at arrival at 
Wind a velocity of 813± 106 km s" 1 (Figure US]) . The average speed of the ICME 
sheath region between the shock and the magnetic cloud as measured in-situ by 
Wind, was « 720 km s^ 1 . Both methods slightly overestimate the impact velocity 
by about 100 km s _1 . In this case, differences in the velocity-time profiles are 
more pronounced: F<f> shows a nearly constant velocity between 800 and 900 
km s _1 while HM delivers an acceleration up to « 100 solar radii to about 1000 
km s _1 followed by a deceleration to 800 km s _1 . Since the irregularity in the 
velocity profile happens in the transition region between the fields of view of HI1 
and HI2 it could also be related to the different sensitivities of the cameras. The 
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Figure 10. Propagation directions for Fixed-</> (top panel) and Harmonic Mean (bottom 
panel). The black crosses show the differences between the calculated distance at arrival time 
and the distance of the in-situ spacecraft from Sun-centre for different propagation angles. 
The red asterisks indicate the same approach but for the difference of the calculated velocity 
and the in-situ measured velocity at arrival time. 



signal of HI1 turns faint rather quickly at larger elongations while HI2 seems to 
overexpose at lower elongations what yields to a discontinuity of the front and 
makes the tracking of the leading edge difficult. These circumstances could be a 
possible reason for the untypical shape of the velocity evolution. 

Table ffl lists the results of the three events for the two conversion methods, 
F</> and HM, and the relevant in-situ measurements. 



5. Discussion and Conclusions 



The purpos e of our stu dy was to introduce a new technique, pointed out by 
Mostl et all (j2009ll20ld) . to additionally use the constraints imposed by in-situ 



measurements to derive ICME kinematics end-to-end from the Sun to 1 AU. 
This method uses single-spacecraft heliosphcric-imager data together with single- 
spacecraft in-situ data and assumes a constant propagation direction. 
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Figure 11. Top: Resulting distance-time profile for F<f>. The crosses are the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal 
line shows the arrival time of the ICME at STEREO-B. Middle: The solid curve shows the 
derivation of the fit, which is done to determine the velocities of this event. The gray area 
indicates the standard deviation of the measurements. The horizontal line indicates the in-situ 
measured velocity at STEREO-B. Bottom: Residuals of the fit and the direct measurements. 



Kinematics covering measurements in the HI1 and HI2 FoV were analyzed 
for three well-observed CMEs on their way from Sun to 1 AU. Images from 
the Heliospheric Imagers onboard the NASA STEREO mission were used to 
produce time-elongation plots ( Jmaps) from which the ICME measurements were 
derived. The measured elongation angles were converted into radial distances 
(in units of solar radii) by using two methods — Fixed-(/) (F^>) and Harmonic 
Mean (HM) — which approximate the ICME front as a point (Ftp) or a circle 
(HM). By combining remote sensing with in-situ measurements it was possible 
to calculate the propagation directions as well as the kinematics of the events 
within the geometrical assumptions of the conversion methods. We constrained 
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Figure 12. Top: Resulting distance-time profile for HM. The crosses are the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal 
line shows the arrival time of the ICME at STEREO-B. Middle: The solid curve shows the 
derivation of the fit, which is done to determine the velocities of this event. The gray area 
indicates the standard deviation of the measurements. The horizontal line indicates the in-situ 
measured velocity at STEREO-B. Bottom: Residuals of the fit and the direct measurements. 



the ICME direction in such a way that the ICME distance-time and velocity- 
time profiles are most consiste nt with in-situ measu rements of the arrival time 
and the velocity on arrival. In lTemmer et al\ (|201li ) the velocity profiles of the 
same three events were compared to th e ambient solar wind, modeled by the 
ENLIL 3D MHD model (|Odstrcill 120031) . 

The ICME of 02-06 June 2008 was clearly identified in the Jmap. It was a 
slow event (~ 390 km s _1 ) that was embedded in the solar wind. It propagated 
toward STEREO-B, where clear signatures of a magnetic flux rope could be iden- 
tified in-situ. The F(j> method reveals a propagation direction toward STEREO-B 
(4>F<p = E24 ±3°), while the HM method assumes a wider longitudinal extent 
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Figure 13. In-situ measured proton density from Wind (top). The first red line from the 
left marks the time of the ICME arriving at Wind and others delimit two strong peaks in 
the proton density. The lower panel shows the Jmap produced from remote sensing data of 
heliospheric images of STEREO-A with overplotted measurement points (red crosses). The 
white horizontal line indicates the position of Wind. 



which may be improbable for this ICME event (4> hm = E51 ± 6°), because the 
exten t in the ecliptic s hould be around 25 degrees ( Wood. Howard, and SockerL 
2010J ). The method of IWood. Howard, and Sockerl (|2010l) where synthetic im- 



ages of a tube-like CME are adjusted to fit the CME's appearance in the 
STEREO/HI imag es, yields a directio n in longitude of E34°, and various di- 
rections derived by iMostl et~gj. I (|2009h range from E24° to E45°, all consistent 
with the directions derived in this paper. 

The event of 13-18 February 2009 was a slow event too (« 320 km s _1 ). Its 
parameters were measured in-situ by STEREO-B but do not show a clear three- 
part structure as for the June 2008 event. This CME originated from the limb 
and seems to have a wider longitudinal extension in the ecliptic and therefore 
was difficult to investigate since it dimmed rather quickly in the Jmap because 
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Figure 14. Propagation directions for Fixcd-</> (top panel) and Harmonic Mean (bottom 
panel). The black crosses show the differences between the calculated distance at arrival time 
and the distance of the in-situ spacecraft from Sun-centre for different propagation angles. 
The red asterisks indicate the same approach but for the difference of the calculated velocity 
and the in-situ measured velocity at arrival time. 



it left the TS early. The results of the used methods differ again by about 25 
de grees (<f>F<t> = E 36±4 °, (pHM =E61 ±3°) and agree with the directions derived 
by iMostl etalj (|201lh who found E35° (F0 fitting) and E63° (HM fitting), 
respectively. The kinematics of both conversion methods are consistent with 
the in-situ measured velocity. 

The event of 03-05 April 2010 is an outstanding event for this solar minimum. 
In contrast to the other two events it was relatively fast (w 835 km s _1 ). The 



derived propagatio n directions r ange f rom ipp^ — W3 ± 4° to 4>hm = E25 ± 10° 
from Earth. Here, Mostl et al. ( 2010l ) found a direction of longitude WO ± 5° 
usin g various m e thods (triangulation, Fixed-0 fitting, and forward modeling), 
and Wood et al. ( 201 lh find a direction of W2°, which makes our Fixed-(/> result 
more consistent with the others quoted in the literature. The velocity profile 
shows a different evolution for both methods. While F(j> results in a nearly 
constant velocity between ps 750 and 850 km s _1 , HM shows an acceleration up 
to « 100 R,0 followed by a slight deceleration. The event seems to be decelerated 
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Figure 15. Top: Resulting distance-time profile for F</>. The crosses arc the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal line 
shows the arrival time of the ICME at Wind. Middle: The solid curve shows the derivation 
of the fit, which is done to determine the velocities of this event. The gray area indicates the 
standard deviation of the measurements. The horizontal line indicates the in-situ measured 
velocity at Wind. Bottom: Residuals of the fit and the direct measurements. 



from the slower ambient solar wind, in which the ICME was embedded (e.g. 



Gopalswamv et all 120011) 



The interpretation of the hcliospheric images is rather difficult. Both conver- 
sion methods (Fixed-;/) and Harmonic Mean) are useful to analyze the kine- 
matics of ICMEs on their way through the heliosphere. How appropriate it 
is to use one or the other clearly depends on the chosen event. The big ad- 
vantage of these methods compared to fitting metho ds such as the Sheeley- 
Rouillard ( Sheelev et al. . 1999t Rouillard et al\. 120081) or the Harmonic Mean 



20irt I.MomI ,1. all . l201lh~ is that a constant velocity 



fitting method (jLueaz . 
is not assumed. On the other hand they are not useable for forecasting be- 
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Figure 16. Top: Resulting distance-time profile for HM. The crosses are the converted mean 
values of the direct measurements and the red line is the spline fit. The blue horizontal line 
shows the arrival time of the ICME at Wind. Middle: The solid curve shows the derivation 
of the fit, which is done to determine the velocities of this event. The gray area indicates the 
standard deviation of the measurements. The horizontal line indicates the in-situ measured 
velocity at Wind. Bottom: Residuals of the fit and the direct measurements. 



cause remote observations over the whole distance and in-situ measurements at 
1 AU are used as input to constrain CME kinematics. T o assess the validity 
of the technique, numerical simulations sho uld be used (e.g. lOdstrcil and Pizzol 
2009; L ugaz. Roussev. and Gombosil 120111) , because there the distance-time and 
velocity-time functions of the ICME front and its shape are known. Another 
disadvantage of the method is that for the calculation of the propagation direc- 
tion only the remotely sensed elongation value at arrival time at 1 AU is used 
and constrained with in-situ data. Assuming constant direction, the resulting 
direction is then used to convert the whole track of the remotely sensed CME. 
Using only one spacecraft from one vantage point it is not possible to derive the 
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Table 1. Summary of the results for all three events. The table lists the 
propagation angle, the arrival time, the mean velocity, and the impact 
velocity of the ICME for the used methods. 



02- 06 June 2008 F0 HM in-situ 

propagation angle E24 ±3° E51 ± 6° 

arrival day 06 Jun. 2008 

arrival time 15:35 

mean velocity [km s" 1 ] 389 ± 48 395 ± 35 

impact velocity [km s" 1 ] 432 ± 61 422 ± 44 400 

12-18 February 2009 F0 HM in-situ 

propagation angle E36 ± 4° E61 ± 3° 

arrival day 18 Feb. 2010 

arrival time 10:00 

mean velocity [km s" 1 ] 310 ± 29 325 ± 28 

impact velocity [km s" 1 ] 416 ± 35 398 ± 34 362 

03- 05 April 2010 F0 HM in-situ 

propagation angle W3 ± 4° E25 ± 10° 

arrival day 05 Apr. 2010 

arrival time 07:58 

mean velocity [km s" 1 ] 829 ± 122 854 ± 100 

impact velocity [km s" 1 ] 825 ± 129 813 ± 106 720 



direction of e very point along the CME-track in contrast to the triangulation 



technique by (jLiu et all l2010a[) 



All this is needed as a basis to be able to better forecast the direction and 
arrival time of coronal mass ejections using empirical or numerical propagation 
models. It will be necessary to investigate a large number of events to reveal their 
different kinematics in order to be able to apply the most appropriate method. 
Furthermore, the geometrical limitations of the Ftp (point shaped CME) and 
HM (circle shaped CME) methods should be adjusted to consider the different 
widths of the CMEs. 
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